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ABSTRACT 

We describe ISO a and BIMA observations of the z = 0.04 Seyfert 1 ring galaxy NGC 985 which suggest 
close parallels with some quasar host galaxies. NGC 985 contains two closely spaced nuclei embedded 
in an R x / 4 -law stellar bulge and an outer ring, evidence of an ongoing merger. The system contains 
~1.8x 10 10 M© of highly disturbed molecular gas which lies in an asymmetric bar-like structure with the 
q . peak in observed CO column densities significantly offset from the compact double nucleus. In contrast 

to this, the ISO observations show strong dust emission centered on the AGN, located in one of the two 
nuclei. Fainter CO, MIR, and radio continuum emission provide a glimpse of the complexities of star 
d • formation in the outer ring. 

An analysis of the kinematics of the main CO emission reveal evidence for two dynamically dis- 
tinct molecular components within NGC 985. The first is a set of isolated super-giant molecular clouds 
(SGMCs) which are concentrated within 9-10 kpc of the active nucleus. Although randomly distributed 
about the center, the clouds may form part of a clumpy highly-disturbed disk which may be either just 
forming around double- nucleus (the making of a ULIRG), or alternatively may be in the process of being 
disrupted-perhaps as a result of a powerful nuclear outflow (the breaking of a ULIRG). A second major 
concentration of CO lies offset from the double-nucleus in a dynamically coherent ridge of emission in 
which powerful star formation is occurring. We tentatively associate CO emission with two out of six 
UV absorption-lines seen in the blue wing of the very broad Lya emission. Such an association would 
imply a complex inter-relationship between the nuclear CO cloud population in colliding systems, and 
AGN-driven winds. 

a ISO was an ESA project with instruments funded by ESA Member States (especially the PI countries: 
France, Germany, the Netherlands and the United Kingdom) with the participation of ISAS and NASA. 

Subject headings: infrared: galaxies — galaxies: individual (NGC 985) — galaxies: interactions — 
galaxies: Seyfert — galaxies: starburst — quasars: absorption lines 



1. INTRODUCTION 

The detection of large quantities of molecular gas in the 
host galaxies of quasars (e.g. Scoville et al. 1993; Schin- 
nerer et al. 1998; Guilloteau et al. 1999; Evans et al. 2001; 
Papadopoulos et al. 2001) lends strong support to the 



idea that quasars and molecular-rich Ultraluminous In- 
frared Galaxies (ULIRGs) may be related in evolutionary 
terms - an idea first suggested on the basis of their sim- 
ilar far-infrared (FIR) properties and comparable space 
densities (Soifer et al. 1986; Sanders et al. 1988; Evans 
et al. 1999). Such ideas find support in the HST imag- 



1 Now at SIRTF Science Center, California Institute of Technology, Mailcode 249-6, 1200 E. California Av., Pasadena, CA 91125 
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ing of nearby quasars, which show that some quasar host 
galaxies have disturbed morphologies, and are sometimes 
"caught in the act" of merging (e.g. Bahcall ct al. 1995, 
1997). Given the potential importance of the molecular 
gas (either for quasar fueling, or for the initiation of major 
star formation activity in the host), we are investigating 
some nearby examples of galaxies which share some molec- 
ular and morphological similarities to host quasars, but are 
close enough that detailed observations can be made. 

NGC 985 (Mrk 1048, VV 285) is a peculiar galaxy which 
was first observed spectroscopically by de Vaucouleurs & 
de Vaucouleurs (1975), and contains a Seyfert 1 nucleus 
in the southern part of the ring. Very broad hydrogen 
and He I A5876 lines in excess of 13,000kms _1 are seen 
in the nucleus, along with other indications of high ex- 
citation suggesting the existence of an optically-thick ac- 
cretion disk (Stanga et al. 1991). The nucleus lies at the 
center of a large IR-bright bulge. Observations at optical 
and near-IR wavelengths have demonstrated that this nu- 
cleus is double (Rodriguez Espinosa & Stanga 1990; Apple- 
ton & Marcum 1993; Perez Garcia & Rodriguez Espinosa 
1996, - see also archival HST WFPC and NICMOS im- 
ages) with a separation of only 3 arcsecs (~2.5 kpc). This 
strongly suggests that the galaxy is undergoing a merger 
in which the nucleus of the intruder galaxy responsible for 
the ring formation is already sinking towards the core of 
the primary galaxy (Perez Garcia & Rodriguez Espinosa 
1996). The one-armed spiral/ring, and the offset nature of 
the nucleus that dominates the luminosity from the radio 
through to the X-rays, suggest some similarities with some 
recently imaged quasar hosts galaxies. In particular, the 
quasars PKS 2349-014 and PKS 0316-346 (Bahcall et al. 
1995, 1997), have very similar morphologies to NGC 985, 
containing highly offset nuclei and outer ring loops. Al- 
though milder examples of offset nuclei are common in 
samples of collisional ring galaxies, offsets as extreme as 
the one in NGC 985 are very rare (see Appleton & Struck- 
Marcell 1996, for a review). 

NGC 985 has an IR luminosity of L IR = 1.8 x 10 11 L , 
which is the largest in a sample of ring galaxies observed 
by IRAS (Appleton & Struck-Marcell 1987), placing it in 
the category of LIRGs (Sanders & Mirabel 1996). The 
absolute magnitude of NGC 985 (My = -22.58 from Ap- 
pleton & Marston 1997) places it in a similar range of 
luminosities to that of radio-loud quasar hosts (Bahcall et 
al. 1997), despite the fact that NGC 985 contains a radio 
source of only a few mJy (this paper - see also Ulves- 
tad & Wilson 1984). NGC 985 is a powerful X-ray source 
(Brandt et al. 1994) with a spectrum which is consistent 
with both a hard power-law component and a substantial 
soft component. The X-ray spectrum is best fitted if an 
additional warm-absorber is placed along the line of sight 
to the AGN. 

Recent HST STIS observations (S. Penton & J. Stocke, 
U. of Colorado: Private Communication) also indicate the 
presence of at least 6 separate narrow UV absorption lines 
in the blue wing of the broad Lya emission line. We will 
argue in this paper that two of these lines may have CO 
emission counterparts, suggesting a possible interaction 
between a nuclear AGN-driven outflow and the molecu- 
lar gas. 



The primary nucleus (hereafer Nl-following the ter- 
minology of Perez Garcia & Rodriguez Espinosa 1996) 
of NGC 985 has a systemic velocity of 12,814±15kms~ 1 
(z = 0.04) based on absorption-line studies performed by 
Arribas et al. (1999). The secondary nucleus (N2) has a 
systemic velocity of 13,096±15kms _1 . Assuming a Hub- 
ble constant of 75 km s -1 Mpc -1 , we adopt throughout this 
paper a distance to NGC 985 of 170 Mpc, so in our figures 
an angular distance of 10 arcsec corresponds to a projected 
distance of ^8 kpc. 

Details on our observations are presented in Section 2 
and the main imaging results are shown in Section 3. In 
Sections 4 and 5 we elaborate on the kinematics and dy- 
namics of the galaxy based on the CO, while in Sections 
6 and 7 we present the MIR evidence for peculiarities in 
the star formation properties of the nucleus. In Section 8 
we discuss the physical conditions in the outer ring H II 
regions, and in Section 9 we explore a potential link with 
of the observed UV absorption lines. In Section 10 we dis- 
cuss the current and future evolution of NGC 985, and our 
conclusions are presented in Section 11. 

2. OBSERVATIONS 
2.1. Space-based Mid-IR Observations 

NGC 985 was observed by ISOCAM, a 32x32 pixel array 
on board the ISO satellite (Cesarsky et al. 1996; Kcssler et 
al. 1996) on December 31 1997 (ISO revolution 776) as part 
of an ISO-GO observation of bright northern ring galaxies 
(see also Charmandaris et al. 1999; Appleton ct al. 1999, 
2000; Charmandaris et al. 2001). The galaxy was imaged 
through 4 filters, LW1 (4.5 [4.0-5.0] /zm), LW7 (9.62 [8.5- 
10.7] /zm), LW8(11.4 [10.7-12.0] /zm), and LW3(15.0 [12.0- 
18.0] /zm) with on-source time of ~10 minutes per filter. 

A lens resulting in a 3 arcsec pixel field of view, was 
used to create a 3x3 raster map in a "microscan" mode 
in both directions (the telescope was moved 4 arcsecs be- 
tween each element of the raster in a regular pattern), 
and the full width at half maximum of the point-spread- 
function (PSF) was ^4 arcsecs at 12 /im. The overall field 
of view was 104x104 arcsecs 2 , which easily encompassed 
the whole galaxy which is over 40 arcsecs across. 

The data were analyzed with the CAM Interactive Anal- 
ysis software (CIA) 2 . A dark model taking into account 
the observing time parameters was subtracted. Cosmic ray 
hits were removed by applying a wavelet transform method 
(Starck ct al. 1997). Corrections of detector memory ef- 
fects were done applying the Fouks-Schubert's method 
(Coulais & Abergel 2000). The sky was subtracted us- 
ing the emission from areas of the map well outside the 
galaxy and its value was within the one expected from the 
coarse DIRBE zodiacal maps. The flat field correction was 
performed using the library of calibration data. Finally, 
individual exposures were combined using shift techniques 
in order to correct the effect of jittering due to the satellite 
motions (amplitude ~1 arcsec). Based on our experience 
with ISOCAM data reduction and application of the above 
methods to a large set of observations, we estimate that 
the absolute uncertainty in our photometry measurements 
is less than 20%. 

2.2. Ground-based Radio and Millimeter Observations 



2 CIA is a joint development by the ESA astrophysics division and the ISOCAM consortium 
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Observations were made with the D-array of the VLA 
on August 25 1992 at 4.9 (A6) and 8.4 GHz (A3. 5 cm). 
These data were flux and phase calibrated, bad data were 
flagged and ignored, and final maps were made using the 
AIPS routine IMAGR. The effects of sidelobes were re- 
moved using the CLEAN algorithm, resulting in images 
of NGC985 with synthesized beam sizes of 19.25x13.23 
arcsecs (A6cm) and 10.67x7.55 arcsecs (A3. 5 cm). 

The 12 CO(1-0) observations of NGC985 were car- 
ried out with the Berkeley-Illinois-Maryland Association 
(BIMA) 10-element millimeter array (Welch et al. 1996) 
at Hat Creek, northern California. The preliminary BIMA 
observations were carried with one full track (7-8 hrs) of C 
array and D array observations in April 18 1998 and June 
22 1998, respectively. A nearby phase calibration source 
0238+166 was observed for several minutes before and af- 
ter every half-hour of observations of NGC985. Further 
observations with both the C array in May 21-24 1999 
(three 6hr tracks and one short 3hr track), and the D 
array in July 15 and 17 1999 (two 3-4 hrs tracks). The 
source 0339-017 was observed as a phase reference instead 
of 0238+166 during this observing because it was closer 
to NGC 985, and was bright enough to be used at this 
time (the source is variable). Individual antenna system- 
temperatures averaged 350 K, and were typically within 
the range 250 to 550 K. Mars and/or Uranus was observed 
for flux calibration. These data were flux and phase cal- 
ibrated, and bad data flagged and channel maps made 
using the MIRIAD software. The synthesized beam of the 
final maps was 11.94x6.81 arcsecs with a position angle 
of -1 degrees. The receivers were tuned to 110.50 GHz 
corresponding to a rcdshift of z = 0.043, and the correla- 
tor was configured to cover a total velocity range of just 
over 1500 km s -1 at a spectral resolution of ~ 8 km s -1 (3.1 
MHz). 

A data-cube containing 40 channel maps was made, 
each channel covering a velocity width and separation 
of 25kms~ 1 , centered on the heliocentric velocity of 
12,934 km s -1 , and so cover a velocity range of almost 
1000 km s -1 . Those maps with bright emission were 
CLEANcd to remove the effects of sidelobes. Subse- 
quent analysis of these data was performed in AIPS and 
MIRIAD. ' 

2.3. Astrometry 

As the next section shows, a large offset is found be- 
tween the peak in the map of CO column densities, and 
the position of the AGN as defined by observations at ra- 
dio and IR wavelengths. Because our observations were 
made under excellent conditions over seven different tracks 
with two arrays, BIMA was expected to return exception- 
ally high-precision astrometry. However to confirm this, 
we constructed maps of the galaxy by splitting the ob- 
servations into two separate independent data-sets which 
used different combinations of phase calibrators and array 
configurations. The results showed identical emission fea- 
tures with excellent positional and flux agreement. We are 
therefore confident that the BIMA astrometry is of very 
high quality. 

Similar arguments can be made for the VLA data which 

3 We also independently checked the astrometry in the 2MASS image 
of the image with an astrometric measurement made in the optical by 



has subarcsec accuracy, and the AGN core was indepen- 
dently observed by Ulvestad & Wilson (1984). Later in 
the paper we also compare our observations with UKIRT 
K-band IR images published previously by Applcton & 
Marcum (1993). These deep observations were registered 
to the astrometrically accurate 2MASS K-band image of 
NGC 985 3 while the details of the registration of the MIR 
images is discussed in the next section. We are therefore 
confident that the registration of all the various datasets 
used in the forthcoming discussion are accurate to an arc- 
second or better - an angular scale much smaller than the 
observed offset in the CO distribution - the result that 
prompted us to perform these checks. 

3. IMAGING RESULTS AT DIFFERENT WAVELENGTHS 

3.1. Visible, IR and Radio Continuum Comparisons 

The images obtained with ISOCAM at wavelengths of 
4.5, 9.62, 11.4 and 15 /im respectively are displayed in Fig- 
ure la-d. For reference, Figure 2a shows a contour map 
of the 15 /im emission superimposed on a grey-scale It- 
band image of the galaxy. We also present in Figure 2b 
the MIR spectral energy distribution (SED) of the nucleus, 
the properties of which are discussed in detail in Section 7. 
The 15 /im filter is the broadest and has the best signal to 
noise ratio of the four filters used. The images obtained in 
the three of the longer wavelength filters (Fig. lb, c, and 
d) show very similar MIR distributions. The strongest 
emission comes from the nucleus of the galaxy and consid- 
erable extended emission is present. Emission is also seen 
near the two optically prominent H II regions in the west- 
ern part of the outer ring. In addition fainter emission is 
detected immediately adjacent to the nucleus to the west, 
especially in Figure lc and d along the optically prominent 
emission ridge. Only the center of the galaxy is detected 
at 4.5 /*m, but observations made with this filter are less 
sensitive than those made at other wavelengths. 

Table 1 lists the MIR fluxes of the four ISOCAM filters 
for the central bulge/nucleus and for the bright "knot" 
seen near the western H II regions in the ring. Additional 
MIR observations of NGC 985 from the ground had also 
been performed by Maiolino et al. (1995) who found that 
the 10.6 /tm (N-band) flux density within a 5.3 arcsec aper- 
ture centered on the nucleus of NGC 985 is 143 mJy. This 
agrees to better than 10% with the average flux density 
we obtained at with the 9.7 and 11.4 /tm ISOCAM filters, 
and is conclusive proof that the majority of the 10 /im flux 
from NGC 985 is very centrally concentrated, since the ISO 
flux was determined over a much larger aperture than the 
ground-based observation and yet the same flux was ob- 
served. This is very different from the situation at 2.2 /im, 
where the AGN's contribution is much smaller than that 
of the stellar bulge (see Section 4) in the inner regions. 

In order to register the MIR maps with the optical, near- 
IR, and radio maps, we could not assume that the absolute 
astrometry of ISO was better than a few arcsecs (unlike 
the relative astrometry which is much more accurate) . We 
have made the reasonable assumption that the bright peak 
in the 15 /tm map corresponds to the Seyfert nucleus de- 
tected at radio wavelengths. We registered this peak to the 

by comparing the position of a foreground star in the northern portion 
(Jeske 1986). They showed perfect agreement. 
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corresponding peak found from the VLA radio maps. Fig- 
ures 3a and b show the radio continuum maps correspond- 
ing to a wavelength of 6 and 3.5 cm respectively. Figure 3c 
shows the 3.5 cm map superimposed on the R-band image 
of the galaxy. The 3.5 cm map, which has the better spatial 
resolution, shows a bright radio source we associated with 
the Seyfert nucleus, and a second component which lies to 
the west just inside the outer ring. The radio fluxes are 
given in Table 1. We fitted a Gaussian to the peaks in both 
radio maps (using the AIPS routine IMFIT) to determine 
an accurate position for the nucleus. Both positions agreed 
within a fraction of an arcsec. The centroid in both maps 
was found to occur at a(J2000) = 02 h 34™37.74! (±0.01!) 
and <5(J2000) = -08?47.'15.0'/ (±0.27'.'). 

One potential difficulty of using the position of the radio 
core for the ISO image registration is the possibility that 
the nucleus might be heavily obscured, even at MIR wave- 
lengths, (e.g. Rigopoulou et al. 1999; Soifer et al. 1999, 
2001), and the peak may not be spatially coincident with 
the radio nucleus. However we can exclude this possibil- 
ity because we find excellent positional agreement between 
the radio nucleus (above) and a K-band image of the nu- 
cleus obtained from the astrometrically accurate 2MASS 
survey (Skrutskie et al. 1997). Since the position of the 
near-IR nucleus and the radio nucleus correspond, we can 
safely assume that obscuration is not a factor in determin- 
ing the position of the AGN core at wavelengths longward 
of 2 fim. 

Figure 4 displays the 3.5 cm radio map superimposed 
on images made at 15 ^m with ISOCAM (Fig. 4a), and 
at 2.2 /tm with UKIRT (Fig. 4b, an image based on near- 
IR data from Appleton & Marcum 1993). This shows 
that the MIR and radio emission have similar distribu- 
tions. Figure 4b shows that the radio source is associ- 
ated with the eastern (brighter) component of the double 
nucleus. No obvious radio emission is found associated 
with the second nucleus, consistent with the suggestion by 
Perez Garcia & Rodriguez Espinosa (1996) that the second 
nucleus is not active 4 . 

3.2. CO Distribution and Integrated Properties of 
NGC985 

Figure 5 shows a map of the integrated 12 CO(1-0) emis- 
sion from NGC985 superimposed on optical (Fig. 5a), 
15 /im MIR (Fig. 5b), and 3.5cm radio maps (Fig. 5c) 
of the galaxy. The emission is distributed in an elongated 
structure with two almost equal peaks separated by 13-15 
arcsecs. A weaker third CO concentration is seen to the 
west, near a part of the ring which contains two powerful 
H II regions. No emission is seen associated with the large 
stellar ring which loops to the north of the nucleus/budge 
at large radii. 

A striking offset of about 5 arcsecs is found between 
the main body of CO emission, and the bright radio/IR 
peak. Indeed the majority of the CO avoids the nucleus 
and is mainly concentrated to the west, where some of it 
is associated with the peculiar ridge of optical/IR emis- 
sion extending away from the center. Such offsets are un- 
usual in collisional galaxies, where either the gas is cen- 

4 P. N. Appleton, in a personal communication, mistakenly reported 
source, a report that was included in their paper on the double nucleus 
radio source appears double, the second component is seen on a much 
The mistake is regretable. 



tered on (or straddles) the nucleus. In the unusual case 
of Arp220, a galaxy with two nuclei probably undergoing 
a major merger, the CO is mainly concentrated between 
the two nuclei in a massive molecular disk (Scoville et al. 
1997), but is still centrally concentrated. 

Very few ring galaxies have been observed in CO with 
high resolution. In a few cases, single dish observations 
have provided some indication of molecular distribution on 
the scale of the large outer rings. The major single-dish 
study of ring galaxies by Horellou et al. (1995) shows that 
some systems (e.g. II Hz 4) have emission concentrated in 
the outer star-forming rings (see also Higdon et al. 2000, 
for AM0644-741). However, only two systems, Arpll8 
and Arp 143, have been published with high spatial reso- 
lution (Gao et al 1997; Higdon et al. 1997). In Arp 118 
most of the CO is found in a powerful star forming ring, 
with little CO(1-0) emission from the center of the galaxy. 
In contrast to this, Arp 143 (Higdon et al. 1997) and the 
Cartwheel (Horellou et al. 1998) have CO concentrated 
towards the center. 

We have constructed an integrated spectrum of the CO 
emission from the galaxy which is shown in Figure 6, along 
with arrows indicating the velocities of the two nuclei. Also 
shown is the range of velocities quoted in the single-dish 
study of NGC985 by Horellou et al. (1995). This under- 
estimates the full range observed, since it shows the range 
at FWHP for some representative spectra (an integrated 
spectrum is not provided in that work) . Their IRAM 30-m 
12 CO(1-0) observations, obtained with a 22 arcsec beam 
and performed over several pointings across the galaxy, 
agree closely in terms of overall velocity spread and line- 
shape with the BIMA spectrum. There is also excellant 
agreement between the BIMA fluxes and the IRAM obser- 
vations. 

Table 2 presents a compilation of the observed and de- 
rived molecular gas properties based on these observations, 
along with published single dish CO and upper limits to 
the H I mass in the system from other work. As the table 
shows, the BIMA observations retrieve 100% of the sin- 
gle dish flux. Adopting the "standard" conversion from 
CO flux to H2 mass (see Young & Scoville (1991)) given 
byM(H 2 ) = 1.18xl0 4 F co D 2 (1+z)" 1 , with D = 170 Mpc, 
then M(H 2 ) is found to be 1.8 xlO 10 M Q . This conver- 
sion factor is highly uncertain, but we adopt a standard 
conversion for comparison with other work. The ratio of 
L m /M(H 2 ) for NGC985 is 10 (see Table 2). This indi- 
cates that the present star formation activity in NGC 985 
is comparable with a local starburst galaxy, but not as 
as high as typical ULIRGs or quasars (Sanders & Mirabel 
1996; Evans et al. 2001). From a star formation point of 
view, NGC 985 has the potential to become significantly 
more luminous and more star formation efficient in the 
future. 

We note that no neutral hydrogen has yet been de- 
tected from NGC 985. An upper limit to the H I mass 
of < 3.6xl0 9 M Q is derived from Heckman et al. (1978) 
after correcting their upper limit for our adopted value of 
the Hubble constant, and assuming that the H I velocity- 

to J.M. Rodriguez Espinosa that the second nucleus was also a radio 
(Perez Garcia & Rodriguez Espinosa 1996). Although it is true that the 
larger scale and is not associated with the second nucleus (Figure 3c). 
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width is similar to that of the CO profile. This leads to 
a M(H2)/M(H I) ratio of greater than 5.5 implying that 
the bulk of the cool interstellar medium (ISM) in NGC985 
is in molecular form 5 . Mirabel & Sanders (1989) have 
shown that M(H2)/(MH I) ratios as large as 4-20 are usu- 
ally found in advanced merger galaxies that have the high- 
est far-infrared to optical flux ratios, are ultraluminous in 
the infrared, and exhibit OH mcgamascr emission. In this 
context, if the upper limit for the HI mass is confirmed, 
NGC 985 could be considered as an ultraluminous infrared 
galaxy in the making. We will discuss this and another 
possibility later in the paper. 

4. CO KINEMATICS - EVIDENCE FOR TWO DYNAMICALLY 
DISTINCT POPULATIONS 

One reason that the main concentration of CO emission 
is not centered on the nucleus of NGC 985 may have to 
do with the peculiar motions present in the CO gas. To 
demonstrate this, we have constructed a data cube which 
represents the observed CO column densities in each of 40 
channels covering a velocity spread of almost 1000 km s -1 
centered on the velocity of the system. 

Figure 7 shows a series of channel maps over the range 
of velocities 12,639kms _1 to 13,089km s 1 in steps of 
25kms~ 1 intervals-the range in which significant CO 
emission is detected. We note that only a subset of the 
full range of velocity coverage is shown here, since the rest 
of the channels show little significant emission above an 
rms noise of 11 mJy beam -1 . The first contour is close to 
the 3a level (30mJybeam _1 ). 

To aid the discussion, we also constructed a composite 
of the various emission features from Figure 7 and super- 
imposed onto the K-band UKIRT contours of the galaxy. 
Because the kinematics of the galaxy is unusually complex, 
we have split the composites into three parts, representing 
the low, medium, and high velocity regimes in the data 
cube (Fig. 8a, b & c respectively). In the following para- 
graphs we will demonstrate that in Fig. 8a, b, the concen- 
trations of CO are somewhat randomly distributed about 
the nucleus showing very little obvious coherence from one 
channel to the next, whereas, Fig. 8c shows emission that 
is much more easily comprehended, as it displays very clear 
channel-to-channel coherence. This is not a signal-to-noise 
effect, since all the features in the maps are significantly 
above the noise, and are identifiable as emission in the 
single-dish observations. Later, we will demonstrate that 
the emission features seen in the first two panels of Fig. 8 
may have distinctly different dynamical properties from 
the more coherent high-velocity emission. 

The lowest velocity channel maps in which CO is de- 
tected are shown in Figure 8a, covering the range 12,764 
to 12,889 kms -1 . At the these velocities, gas is seen cov- 
ering the eastern AGN (Nl), but shifts away, first to the 
south (V =12,789kms _1 ), and then to the west at veloc- 
ities higher than 12,814kms _1 . Indeed at velocities close 
to the systemic velocity of Nl (12,814kms _1 see Arribas 
et al. 1999) the CO emission is offset from Nl by a few arc- 
sees and does not cover either Nl or the western nucleus 
N2. The tendency of the brighter emission to be offset to 

5 Contrary to NGC 985, the Seyfert galaxy NGC 1144 (=Arp 118) Bra 
in that galaxy were highly peculiar as a result of a violent interaction, 
of at least two lower than in this case. 



the west of the AGN is the reason why the peak in the 
integrated map in Figure 5 is offset from Nl. 

In Figure 8b we show the gas at intermediate velocities. 
Here again, as in Fig. 8a, emission is seen which has little 
obvious spatial coherence from one channel to the next. 
In some channels, the gas is roughly associated with the 
ridge of starlight seen by Rodriguez Espinosa & Stanga 
(1990) extending from the nucleus to the west. We note 
that in one channel (V = 12,939 kms -1 ) there is a com- 
ponent of the CO seen to the northwest of the AGN. This 
corresponds to the position of the possible MIR arc dis- 
cussed later in Section 6, although this could be coinciden- 
tal. However, for the most part, the intermediate velocity 
gas is only loosely associated with the linear optical ridge 
which dominates the western side of NGC 985. 

The high velocity emission from NGC 985 begins to 
show more noticable channel-to-channel coherence. This 
is one of the features that distinguishes it from the lower- 
velocity emission. At V = 12,989 kms^ 1 , gas is seen over 
a wide range of positions, extending all the way from the 
western extreme edge of the ring to the nucleus Nl (this 
also shows that the CO emission spectrum in the direc- 
tion of the nucleus contains multiple components). Emis- 
sion over successively higher velocities is seen to shift pro- 
gressively towards the outer ring, bending northwards (at 
V = 13,039kms _1 ) into the region of two powerful H II 
regions. In these channels, a fainter "finger" of emission 
extends back towards the double- nucleus. At extremely 
high velocities the emission is concentrated at the end 
of the stellar ridge. The highest velocity gas, and that 
with the highest velocity dispersion (150 kms -1 ) is seen 
approximately 14 arcsecs west of the nucleus at the po- 
sition a(J2000) = 02 h 34™36.9!, <5(J2000) = -08?47.T8'/. 
This position corresponds to a region in the linear optical 
emission ridge which exhibits the strongest Ha emission 
as discussed by Rodriguez Espinosa & Stanga (1990) and 
Appleton & Marcum (1993). 

The molecular gas distribution and kinematics at first 
seem quite complicated and difficult to understand. How- 
ever, some insight into the behavior of the clouds comes 
when we realize that there may be two different kinds of 
dynamics being exhibited here: a) cloud-complexes that 
are distributed at random positions close to the nucleus 
and have no obvious optical/IR counterpart, and b) cloud 
structures that have well-defined spatial-velocity coher- 
ence. These latter clouds tend to be associated with the 
obvious optical/IR features, and look similar to the struc- 
tures we are familiar with in studying gas in galaxies (i.e. 
spiral arms, rings, tidal loops etc). 

Many of the individual "clouds" we refer to here are 
quite massive, a few x10 8 Mq, and we will refer to them 
as Super Giant Molecular Clouds (SGMCs). To show that 
two dynamically different CO emission features may be 
present, we show in Figure 9 two ways of visualizing them. 
Figure 9a shows the distribution of position angles of the 
SGMCs (north through east) relative to the AGN (nu- 
cleus Nl) as a function of the radial distance from Nl. 
These are the coordinates of the centroids of each SGMC 
measured channel-by-channcl from the channel maps. The 

sford et al. (1999) showed the H I and CO distribution and kinematics 
Even in that case, the observed M(H2)/M(H I) ratio was 2.6 - a factor 
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centroid of each SGMC was determined (to within 1-2 arc- 
sees = 0.8-1.6 kpc) if it was represented by two or more 
contour levels (>4.5cr) in the channel maps of Fig. 7. In 
some cases where extended emission is observed, only the 
brighter obvious clumps were extracted. However, in most 
cases the identification of a SGMC feature was straight- 
forward and unambiguous. From Figure 9 it is obvious 
that within a radial distance of 9 kpc (10.9 arcsecs) of the 
center, the SGMC complexes are distributed in a random 
manner with respect to the nucleus, whereas, the average 
nature of the emission changes and becomes more ordered 
and coherent at larger radii. This transition was noted in 
the discussion of the previous figure, but is put on a more 
quantifiable footing here. 

The separation of the SGMCs into the two components 
is not restricted to their different spatial distributions, but 
also is apparent in their different kinematic behavior. In 
Fig. 9b, we plot the velocity of the SGMCs (relative to 
the AGN systemic velocity of 12,814kms _1 ) as a func- 
tion of radial distance from the nucleus. It can be seen 
that within 9 kpc (representing the area of the "randomly 
distributed" SGMCs) the velocities of the clouds are dis- 
tributed equally on either side of the systemic velocity of 
the active nucleus Nl, and the clouds appear to dynami- 
cally related to the nucleus. However, for the emission at 
larger radii (mainly those associated with the optical/IR 
ridge feature) the clouds cluster together in one region of 
the diagram, a result of spatial and kinematic coherence. 

The implication of Fig. 9 is that the CO emission in 
NGC985 exhibits two distinct kinds of dynamics. A ran- 
domly distributed (radius < 9-10 kpc) SGMC system, cen- 
tered both spatially and in velocity-space on the nucleus 
Nl, and a more distant coherent SGMC system which is 
highly correlated with obvious optical/IR galaxy compo- 
nents. 

4.1. Is the "random component" rotating? 

It is common when analyzing data cubes containing 
spectral line data to create "moment" maps (isovelocity 
and velocity dispersion) which are often useful in helping 
the astronomer interpret their data. Although we con- 
structed such maps early in our analysis of these data, we 
felt the maps were misleading because, i) the data cube 
contains many places where multiple velocity components 
are seen along the line of sight and, ii) the cube contains 
isolated cloud components which are not obviously related 
to another. This stochasticity in the velocity field leads to 
a very confused pattern in the isovelocity map, making 
it difficult to interpret. Despite these reservations, some 
limited inference can be drawn if we divide it up into two 
velocity regimes. 

In Fig. 10a and b we show the isovelocity maps for the 
cube covering the velocity range 12,684 to 12,934 km s _1 , 
and 12,984 to 13,034 km s -1 respectively. This sepa- 
ration roughly corresponds to the CO emission asso- 
ciated with the randomly distributed nuclear-centered 
clouds (Fig. 10a), and higher- velocity more coherent clouds 
(Fig. 10b) discussed earlier 6 . Despite some major excur- 
sions from regularity, Fig. 10a shows that the isovelocity 
contours for the nuclear-centered low-velocity system may 



show the signature of rotation or gas streaming motions, 
even though they are observed as discrete clouds with 
somewhat random positions relative to the nucleus. The 
higher velocity emission is associated with mainly identi- 
fiable optical/IR features in the outer galaxy and can be 
seen in Fig. 10b. Its pattern confirms our earlier observa- 
tion that the gas in the higher velocity channels behaves 
differently from the low-velocity emission - the steepest 
velocity gradient being associated with the outer eastern 
ring. 

5. THE DYNAMICS OF THE LOW-VELOCITY CO CLOUDS 
AND THE IONIZED INNER DISK OF NGC 985 

In the previous section we presented evidence that the 
CO emission seen within 10 kpc of the double nucleus of 
NGC 985 is composed of bright clumps of emission ran- 
domly distributed about the center, but showing a gen- 
eral east to west drift (as is demonstrated in Fig. 10a) 
which might be interpreted as rotation, or at least coher- 
ent streaming motions. What is the dynamical state of 
this cloud system and is it coupled to the motions of the 
ionized gas and stars seen rotating in quite orderly motions 
near the nucleus of NGC 985? 

Arribas et al. (1999) showed that both the stellar and 
emission line velocity field in the central few arcsecs 
around the nucleus of NGC 985 is consistent with an or- 
derly rotating disk , with a kinematic axis along a PA of 
-45 degrees (this is the same positional angle that connects 
the two nuclei). It is already clear from the isovelocity map 
of Fig. 10a that the east-west drift seen in the CO velocity 
field does not share the same kinematic major axis as the 
ionized gas. Indeed, if we ignore the obvious irregularities 
in the CO velocity field, and naively interpret Fig. 10a as 
indicative of a rotating disk, it would have a kinematic ma- 
jor axis of approximately -90 degrees-quite different from 
the kinematic major axis of the more-uniform ionized disk. 

We can explore the possible connections between the CO 
emission and the ionized disk more clearly in Fig.ll, which 
shows the position of the CO emission centroids of eleven 
of the centrally concentrated SGMCs (and their velocities) 
superimposed on the isovelocity "spider" diagram of the 
ionized gas from the work of Arribas et al. (1999). This 
figure shows the east- west drift in the SGMC positions 
with increasing velocity, but confirms that the clouds are 
definitely not co-rotating with the ionized gas disk. For 
example, if the CO clouds were simply clumps of molec- 
ular material moving in circular orbits with the ionized 
gas, then each CO centroid should lie on a corresponding 
isovelocity line as defined by the velocity field of the ion- 
ized gas. In contrast to this, most of the CO centroids 
lie at velocities which are discrepant from the velocity of 
the ionized gas at the same position. A number of the 
CO centroids lie outside the region where the ionized gas 
velocity field was measured, but even in those cases, the 
velocities corresponding to those CO centroids are incon- 
sistent with any simple extension of the "spider" diagram 
of the ionized gas. The CO clouds arc most discrepant 
to the east and south of the nucleus, where they have ve- 
locities which are approximately 100 km/s too low for the 



6 We stress that even this simplistic division of the cube obscures some of the details that are best seen in Fig 8. For example, the isovelocity 
contour at 12,934 kms -1 in Fig. 10a does not smoothly join in a coherent way with the contours at 12,984kms~ 1 shown in Fig. 10b - a result 
of multiple components of emission along several lines of sight. 
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corresponding ionized gas velocity. Only one SGMC (at 
V = 12989 km/s) has the same velocity as the ionized gas 
at that same position (it may be coincidental-but this is 
close to the position of the second (non-active) nucleus). 
The CO emission, while exhibiting the same general trend 
of east-west motion, appears only loosely coupled to the 
ionized disk, and shows some highly discrepant velocities, 
especially to the south and east. 



5.1. An Unstable Collapsing Disk or Clouds Dominated 
by Random Motions? 

In this section we will explore two different dynamical 
states for the inner CO cloud distribution. The first will 
assume that the east-west drift seen in the CO centroids 
is due to rotation, and we will show that these rotational 
velocities are insufficient to stop the CO from sinking to- 
wards the center when compared with the likely mass in 
old stars present in NGC 985. A second alternative also ex- 
ists, that the clouds travel on non-circular orbits through- 
out the center- allowing for the possibility that they are 
marginally stable. The latter, if true, implies that the cen- 
tral CO distribution is highly disturbed, perhaps a result 
of interactions between the CO clouds and an AGN wind- 
an extra complication that must be taken into account 
when modeling the evolution of gas in merging gas-rich 
systems. 

If we were to interpret the overall east-west trend of ve- 
locity with position across the center of NGC 985 in the 
CO clouds as being rotation, we can ask what mass would 
be required to keep the clouds in stable circular motion. 
Taking the CO centroids both east and west of the nu- 
cleus to be representative of the maximum velocity spread 
(275 kms -1 ), this implies that the mass needed keep the 
clouds in circular orbits at a radius of r=10 arcsecs (8.25 
kpc) would be 2.89 x 10 10 (cosec i) 2 M Q , where i is the in- 
clination of the disk. Assuming an average inclination i = 
45 degrees, this would imply a total mass within 8.25 kpc 
needed to produce circular orbits of M c i rc = 6x 10 10 M . 
We will show, in what follows, that this mass is unrealis- 
tically small, implying that if the clouds are in rotation, 
they have insufficient velocity to remain in circular orbits, 
and are likely to be falling inwards towards the center, 
a result that has implications for the future evolution of 
NGC 985. 

Before we estimate the minimum gravitational mass in 
the inner region of NGC 985 using the flux in K-band light, 
we will make another estimate of the mass required for an 
equilibrium configuration for the CO clouds which does 
not make explicit assumptions about the shape of the or- 
bits. Such an approach is justified if we believe (as con- 
cluded above) that the CO clouds may not be in circular 
orbits. The method, called the projected mass estimator, 
was devised by Bahcall & Tremaine (1981) for the prob- 
lem of the mass of dark matter in small groups of galaxies. 
This may be more appropriate for a system in which we 
have only a small number of "test particles" (the SGMCs) 
moving in the potential of the galaxy. The methods finds 
a weighted mass for a system assuming that the "test par- 
ticles" probing the mass are gravitationally bound to the 
central mass. Making the least assumptions about the el- 
lipticity of the orbits, the relevant equation is: 



24 N 

i=l 

where v Z i is the difference in the recessional velocity be- 
tween the i th SGMC, and the central massive object. Ri is 
its projected distance from the center, and N is the total 
number of SGMCs. 

Putting in the values for the randomly distributed 
SGMCs in the above equation, and assuming the velocity 
and position of the AGN are the barycentric coordinates, 
the total projected mass is found to be M o =1.43xl0 11 
Mq. This is more than twice the mass found assuming 
circular orbits, but uses all the cloud velocities rather than 
just those at the velocity extremes of the distribution, as 
in the circular velocity case. 

How do these dynamically estimated masses compare 
with the mass in stars in the inner regions of NGC 985? 
We can determine a lower limit to the mass of stars by 
estimating the mass associated by old stars via the re- 
lationship developed by Thronson & Greenhouse (1984). 
This approach uses the K-band luminosity of a galactic 
region (they applied it to the center of M51) calibrated 
via measurements of the mass in old stars in the solar 
vicinity. Normally such a calculation could not be applied 
to a galaxy dominated by a powerful AGN, but Apple- 
ton & Marcum (1993) showed that at 2.2 /xm, the bright 
nuclear source is only a minor component of the K-band 
light-which is dominated by old stars in the bulge. We are 
therefore able to remove the contribution of the AGN from 
our calculation. We emphasize that this is likely to be an 
underestimate of the true mass, since it neglects younger 
stars and the possible contribution of dark matter. 

Using the maps and photometry of Appleton & Mar- 
cum (1993), we are able to conclude that within a diame- 
ter of 22 arcsecs (18 kpc) the total K-band flux is 33mJy, 
of which only 13mJy (39%) originates within 2x2 arc- 
secs of the AGN. Hence the total K-band flux from the 
region of the random SGMC component is approximately 
20mJy. This corresponds to a total mass in old stars of 
Moidstars = 1.5xlO n M Q . We consider this a conservative 
lower limit to the true mass. 

This lower limit to the stellar mass in the inner regions 
of NGC 985 is larger than the mass estimated from the ve- 
locity and positions of the CO clouds assuming reasonable 
circular orbits, M id s tars/M c i rc > 2.5, and is closer (but still 
larger than) the mass found assuming the clouds travel on 
non-circular orbits. This implies that either the gas is spi- 
raling into the center on a relatively rapid time-scale, or 
that the gas clouds are exhibiting high peculiar motions- 
perhaps some infalling and some expanding outwards from 
the center. In both cases, this argument, when combined 
with the obvious discrepancies with the velocities of the 
regularly rotating ionized gas disk suggests that the CO 
clouds are in a highly disturbed state and are probably 
not in equilibrium with the gravitational potential. 

6. NUCLEAR ASYMMETRIES - FURTHER EVIDENCE FROM 

ISO 

The CO distribution and kinematics are not the only 
features of this galaxy which suggest that the inner galaxy 
is in a non-equilibrium state. In Figure 12a we show that 
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the MIR emission near the nucleus is not symmetrical, but 
that there is a bright arc or shell to the north and west of 
the nucleus at an angular distance of about 5 arcsecs (see 
also Appleton et al. 2000). The figure shows the applica- 
tion of a simple "unsharp mask" technique to the A15/im 
ISOCAM image. This filtering technique is sensitive to 
small-scale structure in the image on the scale of a few 
arcsecs. 

To further explore this feature, which is also obvious 
on independent images taken at other wavelengths (the 
9.62/zm and 11.4/xm images), we performed a series of ex- 
periments using the theoretical monochromatic ISOCAM 
PSF for A15/im, smoothed appropriately to take into ac- 
count the width of the filter. The PSF was also rotated to 
match exactly the micro-scan direction of the observations. 

We show in Figure 12b the result of an iterative re- 
moval of the PSF from the peak of the Seyfert nucleus. 
The image shows the result of three iterations of removal, 
in which the PSF was subtracted from the peak in the map 
with a loop-gain of 0.5 (similar to the technique CLEAN 
of Hogbom 1974, used in radio astronomy) . The following 
conclusions were drawn from our experiments: a) The nu- 
clear source itself is slightly extended on the scale of a few 
arcsecs and, b) there is an excess of emission to the north 
and west of the nucleus which cannot be explained by the 
properties of the detector or the observing technique. Fur- 
thermore, these bright regions fall on the null of the first 
Airy ring of the optical system - again an indication that 
the emission is real. 

In conclusion, we believe that the ISO images show, on 
several independent filters, evidence of an asymmetry in 
the dust distribution near the nucleus which we are unable 
to explain as an instrumental artifact. The emission may 
be shell-like, but because the emission falls so close to the 
bright nucleus we would prefer to await higher resolution 
observations before making a definitive statements about 
its morphology. However, the emission does support the 
view that the dust and gas are highly disturbed and not 
always directly coupled, since this asymmetric dust feature 
is not seen in the same direction as the main CO distribu- 
tion (although, as discussed in Section 4, there may be a 
faint CO counterpart seen in one channel). 

7. MID-IR SPECTRAL ENERGY DISTRIBUTION 

It has been shown (Lutz et al. 1998; Laurent et al. 2000) 
that a galaxy hosting a dominant AGN is clearly differ- 
ent in the MIR from one dominated by a nuclear star- 
burst because the AGN heats the dusty torus surrounding 
it. Perhaps because of grain destruction in the power- 
ful AGN radiation field, the well known emission bands 
at 6.2, 7.7, 8.6, 11.3 and 12.7 ^m (usually attributed to 
Polycyclic Aromatic Hydrocarbons) are extremely weak 
or completely missing from the MIR spectrum of powerful 
AGNs, and this can provide an additional diagnostic if the 
MIR extinction of the nucleus is not too great. 

To examine the MIR characteristics of the nucleus of 
NGC 985, we display in Figure 2b the four MIR flux densi- 
ties measured over an aperture 12 arcsecs in diameter cen- 
tered on the radio continuum peak (from Table 1). This 
aperture covers the majority of ISO emission associated 
with inner regions of NGC 985. The width of the cor- 
responding ISOCAM broad-band filters is also displayed 



with the horizontal lines. For comparison, we also show 
the MIR spectra of two very different emission regions, 
an AGN dominated MIR spectrum from the nucleus of 
NGC 1068 (see Fig. 4 of Le Floc'h et al. 2001), and a typi- 
cal star formation-dominated MIR spectrum from knot B 
of the Antennae galaxies (see Vigroux et al. 1996; Mirabel 
et al. 1998). Both spectra have been rcdshifted to the same 
rest wavelength as NGC 985 and their flux was normalized 
to the 11.4 ^m flux of NGC 985. 

An inspection of the SED in Figure 2b suggests that 
the longer-wavelength spectrum of the nucleus of NGC 985 
looks more like a nucleus dominated by star formation, 
than one where the dust is heated directly by an AGN, 
since the spectrum of NGC 985 does not rise as fast as 
that of NGC 1068, especially in the 10-15 /xm range. Only 
at A4.5/um is there a hint of a "warm bump" from the 
slightly elevated flux at that wavelength. However, since 
this could be due to contribution from old stars of the 
bulge and we lack more definitive diagnostics in the 5- 
10 ^m range (Laurent et al. 2000; Le Floc'h et al. 2001), 
further observations would be necessary to confirm this 
possible result. 

It is of interest to note that our first-order deconvolution 
of the nuclear emission from the AGN (see Section 6) sug- 
gested that the nuclear emission maybe slightly extended 
relative to the theoretical PSF for the telescope. If so, 
we tentatively associate this emission with the stellar and 
hot-ionized gas disk which is seen in the central few arcsecs 
of the galaxy by Arribas et al. (1999) from optical obser- 
vations. This inner rotating "disk" must be composed of 
both stars and ionized gas since both are seen to co-rotate 
along a line joining the two embedded galactic nuclei. The 
MIR observations may be pinpointing the formation of a 
proto-disk at the core of what will eventually become a 
collisional remnant. Could this be direct evidence for the 
formation of disky cores in elliptical galaxies? Higher res- 
olution MIR observations should allow us to determine 
whether the MIR emission docs indeed coexist with the 
inner rotating ionized disk seen by Arribas et al. (1999). 

8. NON-NUCLEAR STAR FORMATION SITES IN NGC 985 

We had seen from the CO, radio, and MIR maps that 
NGC 985 contains a bright ridge of stars, gas and dust 
that extends from the double nucleus all the way out to 
the stellar ring to the west. Furthermore, in the paper 
by Rodriguez Espinosa & Stanga (1990) it was shown 
that the entire length of the optical/IR ridge is populated 
with powerful H II regions, the most powerful of which 
lie towards the end of the ridge, and in two additional 
H II regions in the ring. Although the ISO observations 
show MIR emission associated with the bright Ha emit- 
ting regions, we note that the MIR emission does not al- 
ways mimic the brightness variations from the optical line- 
emitting regions. For example, the two bright H II regions 
in the outer ring (knot A and B in Rodriguez Espinosa 
& Stanga 1990) are much stronger in the ISO observa- 
tions than the end of the linear structure (knots G & H 
of Rodriguez Espinosa & Stanga 1990), despite the fact 
that both sets of regions have very similar observed Ha 
luminosities. Such differences may relate to differences in 
the relative populations of small and large dust grains in 
the various regions (as suggested by Wilson et al. 2000, 
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to explain differences in the MIR properties of the super- 
giant CO complexes found associated with star forming 
regions in the Antennae galaxies), or they could simply 
reflect different internal optical extinction within the op- 
tical nebulae. 

One explanation for the high star formation rates seen 
along the linear ridge is suggested by the high velocity 
dispersion in the CO gas in this region (see discussion in 
Section 10). For example, Stanford et al. (1990); Wil- 
son ct al. (2000); Gao et al (1997) have shown that in 
NGC 4038/4039, one of the most powerful star formation 
sites corresponds to a region of high CO velocity disper- 
sion (the region called SGMC 4-5 by Wilson et al. 2000). 
In our case, we observe a similar phenomenon. In par- 
ticular, the region towards the end of the linear structure 
corresponds to the highest velocity dispersion seen in the 
molecular gas (150 kms -1 ). This region also contains some 
of the brightest star formation sites seen along the opti- 
cal emission ridge, since knots G and H from Rodriguez 
Espinosa & Stanga (1990) have combined luminosities in 
the Ha line of 7.1xl0 40 crgs~ 1 , a value approximately 15 
times the luminosity of the 30 Doradus nebula in the LMC. 

Even the outer ring of NGC 985 has some odd character- 
istics. For example, we note from the overlays in Figure 2a, 
that although the MIR emission at A15/zm falls almost on 
top of the contours of bright ring H II regions A and B 
of Rodriguez Espinosa & Stanga (1990), the CO emission 
(see Figure 5) falls significantly inside the ring, suggesting 
that the CO near the H II regions themselves has been 
destroyed, but the dust grains associated with the H II re- 
gions, clearly seen in the ISO maps, have not! This under- 
lines the complex relationship between CO emission and 
dust emissivity in colliding galaxies - a point already made 
cogently by Wilson et al. (2000) for the NGC 4038/39 sys- 
tem. 

To further complicate the star-formation story for 
NGC 985, we note that the radio continuum emission from 
the ring also fails to correlate spatially with either the opti- 
cal, radio or the CO emission! The main A3. 5 cm emission 
regions, seen in Figure 3, lies inside the bend in the ring, 
and the offset is significant. This offset between the radio 
continuum and the current star formation sites has been 
seen before in ring galaxies (e.g. Applcton ct al. 1999) and 
may be the natural result of the propagation of a star for- 
mation wave through the outer disk of the target galaxy. 
The radio continuum emission, if created by cosmic rays 
from supernova explosions, will tend to lag behind the 
regions of current star formation site because time must 
elapse for the massive stars to explode. 

We can attempt to estimate the number of supernova 
needed to create the observed radio continuum emission. 
Assuming that the the radio emission is predominantly 
non-thermal (the spectral index is -0.6 which suggests mix 
of thermal and non-thermal emission) then wc can esti- 
mate the supernova rate from equation 8 of Condon & 
Yin (1990) : 

WWHz- 1 ) = 1.3 x lO^^-r-^) (2) 

where Lnt is the non-thermal radio luminosity, a is the 
radio spectral index, and Rsn is the rate of type II super- 
nova. For the western radio component, the radio lumi- 



nosity at 6 cm, L(6cm) = 3.3xl0 21 W Hz^ 1 , which corre- 
sponds to a maximum (assuming all non-thermal emission) 
SN rate Rsn = 0.067yr _1 . Is such a supernova rate rea- 
sonable, given the know rate of star formation in the giant 
H II regions near the radio component? 

To test this, we can calculate the rate of type II super- 
nova that would be expected from knots A and B combined 
(the adjacent knots to the radio emission region), based on 
the star formation rates in those H II regions. If we as- 
sume a Salpeter IMF and a reasonable cutoff to the upper 
mass limit on the IMF we reverse equation 12 of Condon 
& Yin (1990) to: 

Rsn = L Hq (W)/(1.1 x lO 36 )^- 1 (3) 

Using the combined values for knots A and B from 
Rodriguez Espinosa & Stanga (1990) of L(Ha) = 9.0xl0 33 
W, we find Rsn=0.01, a factor of six times too low to 
explain the observed radio emission. However, when cor- 
recting for an internal extinction of 2 magnitudes based on 
optical spectra of the knots (Bransford ct al. 1998) these 
rates come into reasonable agreement. We conclude that 
a burst of star formation in the past, with a magnitude 
similar to that currently occuring further out in the ring is 
quite capable of supplying the required cosmic rays needed 
to power the radio source. 

9. UV ABSORBERS ASSOCIATED WITH CO EMITTING 
CLOUDS: COINCIDENCE OR PHYSICAL CONNECTION? 

Several strong absorption features are seen in the blue 
wings of the broad Lya and N V 1238/1242 (S. Pcnton & 
J. Stocke private communication). Intrinsic lines of this 
kind are found in roughly 15% of Seyfert galaxies, and 
seem to represent a population of warm clouds expelled in 
a nuclear wind. Independent evidence for warm absorbers 
in NGC 985 are found from the X-ray analysis of Brandt 
et al. (1994). The HST spectrum kindly shared with us by 
S. Penton and J. Stocke show six absorption lines seen in 
Lya. All except the most blucshifted line is also seen in ab- 
sorption against N V. The velocities for the lines are listed 
in Table 3. The highest velocity line, around 13,065 km s _1 
in Lya, is narrow and isolated (FWHM of approximately 
50kms _1 ) whereas the other lines, although also narrow, 
may be part of a broad blended system extending over 
a full width of 600 km s -1 . Nevertheless, these individual 
narrow features are (except the most blueshifted one) seen 
in absorption against N VA1239,1243 as well. 

Interestingly, we find that the two most red-shifted lines 
correspond to velocities at which CO gas is seen projected 
against the nucleus. In fact, from Fig. 8 it can be seen 
that, because of the off-set in the main CO distribution, 
gas is seen against the nucleus in only three regimes - be- 
tween V=12664-12689, V=12,764-12789 and V=12989- 
13039 km s -1 . The highest and lowest velocity CO emis- 
sion feature seen against the nucleus correspond quite well 
to the velocities of the warm absorbers. 

The high velocity absorption line (and the CO in emis- 
sion) lies close to the systemic velocity of the non-active 
nucleus N2. The narrow width of the UV absorber implies 
that it might be ISM material associated with N2 (despite 
this being seen in absorption against Nl). 

The UV absorber system observed at V= 12,670 km s -1 
also seems to correspond to the velocity of gas seen in 
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CO emission in front of the nucleus - in fact it cor- 
responds to CO seen in two adjacent channels (spread 
over 50km These emission features form part of the 
randomly-distributed cloud population discussed earlier. 
Could the UV absorber and the CO emission be related 
- perhaps through an interaction of the AGN wind with 
the infalling clouds? One argument against this is that the 
CO clouds are cold, and yet the UV absorption-line clouds 
are warm. However, if an infalling cloud was, through the 
chaotic nature of the infall, to encounter an AGN-driven 
wind, it may become heated and its outer envelope may 
well become ionized. Although our CO observations cover 
two more of the uv-absorption line velocities, no emission 
was detected from the more blue-shifted components. Per- 
haps dense molecular clouds encountering the wind do not 
survive long, and are evaporated rapidly as they become 
entrained in the wind. We note that this may be the first 
time there had been a possible identification between in- 
trinsic Lya absorption line components and CO emission, 
although CO absorption lines have been discovered associ- 
ated with intervening galaxies in the direction of millime- 
ter bright high-z sources (e.g. Wiklind & Combes 1994). 

10. NGC985: PRESENT AND FUTURE PROSPECTS 

The presence of a faint extended (r = lOkpc) and prob- 
ably unstable population of SGMCs in the inner regions 
of NGC 985 is generally unusual for merging or collisional 
systems, although it does share some similarities with 
some well studied cases. For example, in the case of 
Arp 220 (Scoville et al. 1997), the CO, although highly con- 
centrated towards the center, involves two embedded nu- 
clei revolving within a (CO rich) rotating disk. Similarly, 
the peculiar ring galaxy Arp 143 (Higdon et al. 1997), like 
NGC 985, also contains two molecular components, one 
concentrated in the ring, and another concentrated in the 
nucleus. Both these systems have highly centralized CO 
distributions in rapid rotation and may suggest one possi- 
ble direction in which NGC 985 is evolving. Alternatively, 
as we shall discuss below, NGC 985 may actually have al- 
ready been through the dense centrally-concentrated stage 
(ULIRG stage-Sanders & Mirabel (1996)), and its diffuse 
clumpy CO distribution may represent the early stages of 
the disruption of the central disk by a combination of AGN 
and star-formation-driven winds from the center. Both 
possibilities could explain the present state of NGC 985. 

10.1. NGC 985: The Formation or Destruction of a 

ULIRG? 

Despite an incomplete knowledge of the details of the in- 
teraction, the observations presented here provide a com- 
pelling case for a bright future for NGC 985! The peculiar 
CO distribution, decoupled as it is from the inner rapidly 
rotating ionized and stellar disk, suggests the gas is highly 
disturbed, and some of it may be either falling towards 
the center (see Section 4), or may be in marginally stable, 
but non-circular orbits. There are however many questions 
that are raised by the observations. What will be the fate 
of the 2x 10 10 M Q of molecular gas over time? What is the 
significance of the intrinsic warm absorbers seen in the UV 
and X-ray spectra of the Seyfert nucleus and their relation- 
ship to the random component of the CO clouds? Finally, 
is NGC 985 likely to become much more IR-luminous in 



the future, evolving towards a ULIRG state, or is there 
evidence that it has already passed through that state and 
that its activity may well be in decline? 

The current observations provide only partial answers 
to these questions. Firstly, in the absence of a disturb- 
ing nuclear wind, the CO reservoir cannot remain out of 
equilibrium for more than a few crossing times, given the 
dissipative nature of gas. Taking the overall scale of the 
CO structure (about 20-25 arcsecs = 16-20 kpc) and di- 
viding by the total velocity-spread in the CO (Fig. 6) of 
400 km s -1 , provides a rough timescale for the dissipative 
collapse of the gas of 45-50 Myrs. If the gas is indeed in- 
falling and we conservatively assume that 2 such periods 
would be needed to allow the gas to settle, then we might 
expect the large-scale accretion of the gas into the center 
to occur over the next 100 Myrs. If this was a smooth 
process, this would provide an average accretion rate of 
200 M Q yr _1 . Given the peculiar distribution of the gas, 
it is unlikely to be a smooth event. If each of the ran- 
domly distributed SGMCs falls into the center separately 
they will each induce a period of sporadic activity as the 
gas is assimilated into the central disk. It is hard to be- 
lieve that under these condition, the nuclear activity in 
NGC 985 will not increase significantly - perhaps further 
increasing the already luminous nature of the AGN. 

What might be the effect of this infall on an already 
active nucleus? One possibility is that some infalling gas 
will become entrained in a wind presently blowing from the 
Seyfert core. There is evidence from the X-ray spectrum 
of a significant soft X-ray component to NGC 985. Some 
of this may relate to the nuclear star formation which we 
believe must be present in NGC 985 from the ISO spectral 
energy distribution, but some may be due to an outflow- 
ing wind. It is possible that CO clouds raining towards the 
center become entrained and accelerated in the outflowing 
wind and as a result become heated-creating warm outer 
atmospheres which we tentatively associate with the UV 
absorption lines and X-ray absorbers. Indeed it would 
seem natural that inflowing material in tidal interactions 
involving active nuclei would inevitably encounter winds 
from the AGN as some stage. This may explain the fact 
that intrinsic UV absorbers are seen only in a sub-set of 
cases in UV spectra of Seyfert galaxies (the major accre- 
tion stage?). Further observations would be needed to see 
if the appearance of UV absorbers in AGN spectra is cor- 
related with the appearance of highly asymmetric CO dis- 
tributions, like that of NGC 985. In any event, our present 
observations provide for the possibility that the tidal fu- 
eling of AGNs must take into account the possibility of 
potentially complex interactions between the infalling gas 
and the turn-on of an AGN wind. 

The possibility that a nuclear wind could interact, or 
even disrupt the formation of a central nuclear molec- 
ular disk leads to another possible path for the evolu- 
tion of NGC 985. Could the SGMCs seen in the inner 
region represent the ending of the ULIRG stage rather 
than the beginning of it? We have mentioned that many 
interacting systems show highly concentrated CO disks, 
and NGC 985 is quite exceptional in that its CO distribu- 
tion is highly asymmetric-indecd the brighter CO is seen 
away from the nucleus in the outer disk region, and the 
CO in the clouds with highly peculiar motions are sig- 
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nificantly fainter. Could these clouds represent the dis- 
rupted remnants of a much more centrally concentrated 
CO distribution -the stage that is assumed to be associ- 
ated with ULIRGs? 

The question of whether NGC 985 about to become a 
ULIRG, or alternatively has recently passed through that 
stage can be tested with future observations. Observations 
of the elliptical-like bulge component which surrounds the 
double-nucleus may help to determine whether the galaxy 
has recently experienced a major episode of global dis- 
tributed star formation-a symptom of an ongoing or post- 
merger history (Schweizer 1995). Also, if NGC 985 is 
emerging, rather than entering a ULIRG stage, it would 
be expected to show a very extensive soft X-ray emission 
halo resulting from a recent violent star formation history. 
High-resolution multi-wavelength IR observations of the 
bulge may also help to determine what contribution the 
bulge makes to the strong extended mid-IR emission seen 
associated with NGC 985. Such observations would help 
to determine whether the R 1 / 4 -law bulge itself contributes 
a significant amount to the total MIR energy budget of the 
NGC 985 system. 

11. CONCLUSIONS 

Observations of NGC 985 with ISO, the VLA and BIMA 
have led to the following conclusions about this system of 
galaxies: 

1) A MIR luminous AGN is centered on the brighter 
component of a double nucleus. The MIR nucleus, which is 
detected in all 4 ISOCAM bands (from 4.5 itm to 15 /zm) is 
consistent with emission from a rapidly rising thermal con- 
tinuum. The slope of the continuum is consistent with dust 
heated predominantly by star formation, rather than the 
AGN itself. We tentatively associate this slightly-extended 
warm nuclear dust emission with the rapidly rotating stel- 
lar and gaseous disk seen at optical wavelengths by Arribas 
et al. (1999). This disk, which rotates along the line join- 
ing the two embedded galactic nuclei, may be forming as 
a consequence of the merger, resulting from inflowing gas 
to the nucleus. 

2) Our BIMA observations of NGC 985 support the 
view that it is in a violent state of merger. The CO 
emission, though covering the nucleus, is distributed in 
a highly asymmetric distribution, with the peak line-of- 
sight-column significantly offset from the compact double- 
nucleus. This distribution resembles the molecular dis- 
tribution of a number of quasars recently mapped in CO, 
many of which have asymmetric and irregular CO distribu- 
tions (Evans et al. 2001). Most of the cool gas in NGC 985 
is contained within this large molecular reservoir, which, 
when adopting standard conversion factors from CO to H2 
masses, contains 1.8xl0 10 M Q . 

3) Although the highest CO column densities are seen 
away from the nucleus, we isolate a population of SGMCs 
which may be mildly collapsing towards the center on a 
timcscalc of 60-100 Myrs. These clouds, which extend 
within 9-10 kpc of the AGN, were identified through their 
random distribution and unusual motions which are dif- 
ferent from the more coherent motions of material fur- 
ther out in the CO ridge. This two-component nature 
of NGC 985, may indicate a transition time in a merging 
system in which the motions of the clouds begin to loose 



their original ordered motions in the collision, and begin 
to rain down on the center through dissipation. Sporadic 
accretion rates are likely to exceed 100M Q yr~ 1 if all the 
molecular material present fall to the center in a few dy- 
namical crossing times. 

4) NGC 985 currently has a low ratio of L m /M(H 2 ) = 
10. This value lies in the range found for low-luminosity 
starburst galaxies, despite its abnormal CO distribution 
and large, as yet mainly untapped, molecular reservoir. 
This low value may indicate that either i) NGC 985 is 
about to enter the ULIRG stage as the huge molecular 
reservoir accumulates at the center, or ii) that the galaxy 
is emerging from the ULIRG phase, and that the action 
of a powerful AGN or nuclear starburst wind is inhibiting 
star formation in the inner regions. 

5) We tentatively associate two CO emission features in 
the direction of the AGN with two out of six narrow UV 
absorption-lines seen in an HST study of NGC 985. One 
of these CO emission features is associated with the pop- 
ulation of clouds that might be mildly infalling towards 
the center, and this raises the possibility that the infalling 
clouds may sometimes become entrained in an outflowing 
AGN wind. How such interactions may affect the fueling 
of an AGN over time is not clear. However, these interac- 
tions may play an important role in the disruption of cen- 
trally concentrated molecular disks in starburst systems 
that contain powerful AGN. 

6) The CO furthest away from the center shows more or- 
dered motions and is associated with the optical/IR ridge- 
probably a spiral arm or ring-segment containing com- 
pressed gas. The tip of the CO ridge, the most distant 
from the double- nuclei, shows the highest velocity disper- 
sion - a property that seems to correlate with a high rate 
of local star formation. 

7) The ISO observations also reveal warm dust asso- 
ciated with the western optical/IR ridge extending from 
NGC 985. The dust is consistent with grains heated by 
bursts of massive star formation which are observed along 
the ridge and into the outer optical ring of NGC 985. As in 
the case of the Antennae galaxies, the MIR dust emission 
does not always scale linearly with other observed prop- 
erties of the star formation regions (e.g. Ha luminosity, 
molecular cloud mass). We find that in one region, MIR 
emission is seen associated with a bright optical H II com- 
plex, but no 12 CO (1-0) emission is observed at that po- 
sition, indicating a complicated relationship between dust 
grain lifetimes and molecular gas lifetimes in massive star 
formation regions. Our radio continuum observations also 
provide evidence of a previous episode of star formation 
which occurred in the recent past inside the ring. 
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Table 1 

Radio and Mid-Infrared Properties of NGC 985 





ISO-LW1 


ISO-LW7 


ISO-LW8 


ISO-LW3 


VLA 


VLA 


Region 


S(4.5/zm) 


S(9.7^m) 


S(11.4/zm) 


S(15/mi) 


S(3.5cm) 


S(6cm) 




mJy 


mJy 


mJy 


mJy 


mJy 


mJy 


Entire Galaxy 


55.2 


113.5 


154.4 


168.9 


2.28 


4.47 


Nucleus/Bulge a 


46.5 


105.6 


142.5 


163.6 


1.49 


3.38 b 


Eastern SF Knot/rinj 


f 0.8 


5.4 


6.2 


3.5 


0.79 


1.09 


a Flux measured in a 


circular apert 


ure of radius 12 arcsecs 


centered on 


the nucleus. 





b Ulvestad & Wilson (1984) quote a 6 cm flux for the nucleus of 2.5 mJy with the VLA in A 
configuration, but mention that this is probably an underestimate of the flux since they note that 
the flux they detected varied when the array was tapered differently. Our observations, which 
were made with many short spacings in the D configuration, confirm that the A-array observations 
missed a substantial fraction of the flux. 

a Thc ISO fluxes were evaluated over the bright MIR emission knot centered at a(J2000) = 2h 
34m 36.5s, <5(J2000) = -8 d 47' 10", whereas the radio emission was evaluated over the larger area 
defined by the radio maps of the second western source seen in Fig. 3. 
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Table 2 

Global properties of NGC 985 



Notes 



Assumed Distance (Mpc) 


170 




m v (mag) 


13.57 


from Appleton & Marston (1997). 


Lir (L ) 


1.8xlO n 


definition of Sanders & Mirabel (1996). 


M(H I) (M ) 


<3.6xl0 9 


from Heckman et al. (1978) a 


S CO =/ Sdv (Jykms- 1 ) 


54 


from Horellou et al. (1995), filled aperture 


Sco=J Sdv (Jykms -1 ) 


54.8 


(BIMA) this work. 


CO Line Width (FWHM, kms -1 ) 


408 ± 15 


this work. 


M(H 2 ) (M ) 


1.8xl0 10 


this work! 3 


Lir/M(H 2 ) (L /M ) 


10 




M(H 2 )/M(H I) 


>5.5 





a Based on Heckman et al. (1978) after correcting for our adopted value of H . 
b UsingaN(H 2 )/Ico = 3xl0 20 mol cm -2 (Kkms -1 ) -1 
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Table 3 

uv (hst/stis) observations of nuclear absorption lines a 



hya h 


Inferred heliocentric cloud velocity 


A 


krns" 1 


1264.787 ± 0.02 


12,113 ± 4 


1265.317 ± 0.02 


12,243 ± 4 


1265.927 ± 0.06 


12,394 ± 4 


1266.396 ± 0.05 


12,509 ± 4 


1267.047 ± 0.02 


12,670 ± 4 


1268.611 ± 0.00 


13,053 ± 4 



a Kindly provided by S. Penton & J. Stocke. 
b The rest wavelength of Lya is 1215. 67 A. 
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Fig. 1. — a) The ISOCAM 5 (im image of NGC985, with an overlay of nine of its own contours spanned equally logarithmically between 
0.04 and 1.56 mjy pixel - 1 . The coordinates are in J2000. b-c-d) Same as for a) but for the 9.62 /im, 11.4 /xm, and 15 (im image of the galaxy 
respectively. The spacing of the nine contours of each image is set as in a) but with limits [0.02 — 2.47] mjy pixel -1 for the 9.62/xm, [0.02 - 
3.1] mjypixel for the 11.4/xm and [0.02 — 3.4] mjypixel - 1 for the 15/xm respectively. 
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Fig. 2. — a) A contour map of the 15 /an image shown in Fig. Id) superimposed on an optical R-band image of NGC 985 kindly provided 
by A. M. Perez Garcia and J. M. Rodriguez Espinosa (IAC, Spain). The positions of the two nuclei Nl and N2 are marked, b) The nuclear 
spectrum of NGC 985 as measured with a circular aperture 12 arcsec in radius. The width of each broad-band filter used is marked with 
a horizontal line. For comparison we display the spectrum of NGC 1068 (from Fig. 3 of Le Floc'h et al. 2001) in solid line, and that of the 
star forming knot B of the Antennae galaxies (see Vigroux et al. 1996) in dashed line. Both spectra have been red-shifted to the distance of 
NGC 985 and normalized to the flux density of NGC 985 ISOCAM LW8 (11.4/Lim) filter. 
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Fig. 3. — a) The VLA 6 cm (C-band) radio continuum map of NGC 985, with an overlay of its own contours. The contours are 0.3, 0.5, 
0.8, 1.2, 1.9, and 3.0 mjybeam - . b) Same in a) but for the 3.5cm (X-band) radio continuum. The contours are 0.2, 0.3, 0.5, 0.8, 1.3, and 
2.0 mjybeam -1 . c) Same contours as shown in b) but overlayed on the R-band image of NGC 985 shown in Fig. 2a. 




Fig. 4. — a) The 3.5 cm (X-band) contours of NGC 985 overlayed on the 15 (im image of the galaxy. The contour spacing is the same as 
in Fig 3b. b) Same overlay contours as in a) but this time on the UKIRT K-band image of NGC 985 (see Figures 1 and 2 in Appleton & 
Marcum 1993). 
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Fig. 5. — a) The integrated 12 CO(1-0) emission in contours, over the optical R-band image of NGC985. The contour levels are 4.0, 5.2, 
6.8, 8.8, 11.5 and 15.0 Jy beam - 1 kms . b) Same as in a) but the background is the 15 fim image of the galaxy, c) Same as in a) but the 
background is the 3.5 cm (X-band) radio continuum map. 
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Fig. 6. — Integrated spectrum of the 12 CO(1-0) emission from our BIMA data. The optical velocities of the two nuclei, Nl and N2, as well 
as the velocity ranges seen by the single dish observations of Horellou et al. (1995) are marked. 
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Fig. 7. — Channel maps of the 12 CO(1-0) emission from NGC985. The velocity of each channel is marked on the upper left corner. The 
contour levels are nxO. 015 Jy beam - where n=-2,2,3,4,5, and 6. Negative contours are marked with dashed lines. The first positive contour 
is close to the 3<r noise level. 
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Fig. 8. — A composite of the various emission features superimposed on the K-band UKIRT image of the galaxy. Because the kinematics 
of the galaxy is unusually complex we have split the figure into three parts, representing the low velocity, medium velocity and high velocity 
regimes in the data cube (a, b & c respectively). The scale in the upper left corner indicates the angular extend of all the features displayed 
(see Section 4). 
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Fig. 9. — a) The position angle as a function of distance from the nucleus Nl of the scattered and coherent SGMCs observed within 15kpc. 
b) Same as in a) but plotting the velocity of the SGMCs with respect to Nl as a function of distance. Note how at larger radii the clouds 
cluster together in one region of each diagram, a result of spatial and kinematic coherence (see discussion in Section 4). 




Fig. 10. — a) The isovelocity field of the gas for the low-velocity system superimposed on a grey-scale representation of the CO integrated 
column density map-helioccntric velocities arc in kms -1 . b) same as a but for the high velocity system (see text in Section 4.1). 
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Fig. 11. — Emission centroid positions and velocities for eleven SGMCs detected in CO, superimposed on the isovelocity map of the ionized 
gas disk mapped via optical emission lines by Arribas et al. (1999). Note the east-west drift of the CO centroids with increasing velocity and 
the miss-match with the velocity field of the ionized gas disk (see Section 5). 
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Fig. 12. — a) The 15 /im image of NGC 985 after the unsharp mask technique has been applied. The asymmetric emission and a bright 
"ridge" to the north and west of the nucleus, b) 15 /an image after a PSF fit and subtraction of the Seyfert nucleus. The field of view is the 
same as in a). Even though less well defined than the unsharp image one can easily see the same extend in MIR emission to the north of the 
nucleus. 



